Topological Symmetry And Existence of Partial Synchronization 
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We study the relationship between the partial synchronous (PaS) state and the coupling structure 
in general dynamical systems. By the exact proof, we find the sufficient and necessary condition of 
the existence of PaS state for the coupling structure. Our result shows that the symmetry of the 
coupling structure is not the equivalent condition which is supposed before but only the sufficient 
condition. Furthermore, for the existence of the PaS state, the general structure is the equal-degree 
random. 
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An interesting form of dynamical behavior occurs in 
the coupled systems when only some of the subsystems 
behave in the same way, which is Partial Synchronization 
(PaS). The behavior has attracted much attention in the 
physical [l[, biological @, ecological and other sys- 
tems. As the coupling strength between the subsystems 
is small enough, the state of coupled system is turbu- 
lence. Increasing the coupling strength, the PaS state is 
often observed before the Global Complete Synchroniza- 
tion (CS). The existence of various PaS states leads to 
the remarkably complex dynamical behaviors 0] in the 
coupled systems. 

Considering the dynamics Q 




X n+ i =F(X n )+er®CF(X n ), 



(1) 



where X = (x 1 , x 2 , • • • , x N ) represents the status of the 
system. The independent state, which is defined as x 1 ^ 
x 2 7^ • • • =^ x , and the CS state, which is defined as 
x 1 = x 2 = • • • = x" , exist for common coupling ways 
@ between the subsystems. However it is possible for a 
given coupled system without any PaS solutions, Fig. Q] 
just shows an example. Substituting all possible 365 PaS 
solutions (e.g., x 1 = x 2 and x 3 ^ x 4 ^ x 5 ^ x 6 ) into Eq. 
([1]) at the coupling structure shown in Fig. [1] will give 
false statements. 

In recent years, it is believed there are some tight rela- 
tions between the topological symmetry of the coupling 
structure and the PaS state. For example, the asymmet- 
ric PaS pattern that disagrees with a symmetrical struc- 
ture has never been observed 0]; the symmetry group 
theory can be used to describe the partial periodical syn- 
chronous state in some regular structures with the same 
symmetry @; all of PaS states corresponding to each 
topological symmetry in a ring were observed etc. 
So one could suppose that the symmetry is the necessary 
and sufficient condition for the existence of the PaS state. 
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FIG. 1: A topological structure of the coupled system without 
any partial synchronous solution. 



For small amount of subsystems, it seems that the 
above relationship comes into existence. The coupled 
structure Fig. [2)Ja) can be represented mathematically 
by the adjacent matrix 



A A = 



Here, node 1 and 2 is symmetric in the structure, and A4 
will be invariable under a permutation transformation 
1 -> 2 or 2 -> 1. The curves di. 2 {e), d 2 , 3 (e) OH are 
shown in Fig. [U (b), e is the coupling strength. The 
synchronous solution Xi = x 2 7^ X3 ^ X4 is observed in 
the region [0.3, 0.45] U [0.7, 1] of e. Thus, the PaS state 
will be achieved with the corresponding symmetry in A± 
among the symmetrical nodes. 

A more complex case is shown in Fig. [3] (a) with 
the same dynamics as above. There are two clus- 
ters and their nodes are denoted by 1,2, ■■■ ,n\ and 
n\ + 1, ni + 2, • • • ,n% + n 2 . Node i (1 < i < n\) is 
coupled to the nodes i — k, i — k + 1, • • • , i + k, ni + i — 
I, n% + i— Z + l, • • • , rii + i+l, and node n\ +i (1 < i < ni) 
is coupled to the nodes n\ + i — k, n\ + i — k + 1, • • • , 
ni + i + k, i — I, i — I + 1, • • • , i + I. Obviously, there is 
"rotate" symmetry in every cluster. The adjacent matrix 
will be invariant under a "rotate" permutation transfor- 
mation in each of them: e.g., it is 1 — > 2, 2 — > 3, • • • , 
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FIG. 2: (a) A simple, but typical, network with symmetries, 
(b) The average distance di,2 and ^2,3 versus the coupling 
strength e. 
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FIG. 3: (a) The scheme of a topological structure with the 
"rotate" symmetry, (b) The variation of the two clusters 01, 
o 2 and the whole system <r as functions of e with ri\ = = 
100, ft = 40, I = 10. 



m — 1 — ► ni, ni — > 1. Wc define a M x M matrix i?jv/, 
where (R m )m,i = 1, (-Rm)m+i = 1 = 1 ' ' 
and other elements of i?M are 0. Thus the permutation 
matrix [l2j of this transformation will be T<j = ©i? n2 , 
where "©" is direct sum of two matrices. Fig. [3] (b) 
shows the time average of the variation of all subsystems 



a(e) = limjv'^oo w En=i V^=i( x ™ ~ ^ Si=i x «) 2 ' 
and of the two clusters <J\{e), 02(e) with n\ = n% = 100, 
k = 40, Z = 10. For e e [0.45, 1], o\ = a% = 0, namely 
Xi = • • • = x ni , x„ 1+ i = • • • = xat. The PaS solution of 
the "rotate" symmetrical nodes is observed. 

In this Letter, we investigate in detail the relationships 
between the PaS solution and the coupling structure. 
The PaS solution is defined as follow: for a dynamical 
system with phase space R"" 1 , a i^-cluster synchronous 
solution is a if-dimensional subspace, denoted by V , of 
i" 1 ™. It can be represented by 



(2) 



X J 1 


■ 2 

= X l i = • 


• = x'i 








x'-' 


= X*2 = • 


,-"2 

• = X'2 
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= X l K = ■ 


■ = X K 



where x l » denotes the 6th subsystem in the ath clus- 
ter and {rii}fL 1 is the size of each cluster that satisfied 



J2? =1 m = N (N > K > 1). The CS and independent so- 
lution are the particular cases of the definition for K = 1 
and K = N. The relationship between PaS solution and 
the coupling structure could be described by two ques- 
tions as follows: 

Question A: If one finds symmetry in a coupling struc- 
ture, can a corresponding PaS solution be obtained? 

Considering the matrix form of Eq. ^ , 



rx = x, vx e v, 



(3) 



where T is a permutation matrix. X G V is invariant un- 
der the transformation T, so V is the invariant subspace 
of T and the eig-subspace of T corresponding to eigvalue 
1. That V is the invariant subspace of the dynamical 
system Eq. ([T]) requires 



CX e v, VX e V. 



(4) 



Next, if there is a symmetry T in structure C, then C 
will be invariant under a permutation transformation T. 
The mathematical representation is 

T- 1 CT = C, (5) 

Therefore, the matrices C and T are commute, or 



CT = TC. 



(6) 



Question A can be represented by a mathematical 
statement as follow: 

If TX = X, 

then Eq. © Eq. (gj 

Multiplying the two sides of Eq. ^ by X € V, we have 

TCX = CTX,VXe V; (7) 

Combining Eq. (J3]) with Eq. ([7]) gives 

TCX = CX,VXe V. (8) 

So, CX is also the eigvector of T with eigvalue 1. And 
then, Eq. (g]) will be satisfied for all X <E V. Thus, we 
conclude that, for a symmetrical structure, the dynamical 
system has a corresponding PaS solution. 

Question B: If one finds a PaS solution, can the cor- 
responding symmetry in the coupling structure be ob- 
served. 

Here, an interesting example is shown in Fig. g] (a). 
Considering two clusters, each one has n subsystems, ev- 
ery subsystem is randomly connected to [pn] + 1 sub- 
systems in the same cluster (p is a probability and [pn] 
means the integer part of pn) and \p r n] + 1 subsystems in 
another cluster (p r is also a probability) . Each subsystem 
in each cluster has equal degree and the connections be- 
tween two clusters are also equal-degree. Fig. |4](b) shows 
the variance of the two clusters (01, 02) and the whole 
system (a) as functions of the coupling strength e in the 
parameters n = 100, p = 1, p r = 0.5. As e € [0.34,0.72], 
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FIG. 4: (a) The scheme of an equal degree random struc- 
ture, (b) The variance of the two clusters ai, 02 and the 
whole network a as functions of the coupling strength e in 
the parameter m = n.2 = 100, p = 0.5, p r = 1. 

a% =02 = and a > 0, the PaS solution is observed. 
Due to the random connections between the subsystems, 
there is not symmetry in the structure. 

Question B can be represented by a mathematical 
statement as follow: 

If TX = X, 

then Eq. (J) =4> Eq. © 

The following relations could be derived 

C*X G V,TX = X,VX G V 
<=^> TCX = CX 

^> TCX = CTX 

^ (r- 1 cr-c)x = o 

or (CT - TC)X = 0. (9) 

Obviously, Eq. © is not equivalent to Eq. ([5]). We can 
conclude that it is possible for a PaS solution without 
any symmetry in a dynamical system. Fig. [4] just gives 
an example. Another important feature of Eq. © is 
that, in fact, the sufficient and necessary condition can 
be drawn from Eq. 

The component form of Eq. ([9]) is 

N N 

C m „x n - do* = 0,VX = (x\x 2 , • • • , x N ) G V. 

(10) 

for Ti m = 1- Relabelling the subsystems in order to col- 
lect together the subsystems in the same cluster, Eq. @ 
will be rewritten in the form 

Cluster! : 
x 1 = x 2 = • • • = x" 1 = y 1 



Clusters : 

X EC! Uk+1 = x^fc=i " fc+2 = ■ • • = x^fc=i = y s 



ClusterK : 

(11) 



And the general form [l3[ of T will be 

K 

T = ®R ni ,i = l,2,--- ,K. (12) 

»=i 

Then the general form of X will be X = 

(y 1 li,r i i,y 2 li,n 2 ) ' ' ' ,y K li,n K ) T where 1 m ,n is an M x 
iV matrix in which every element is 1. Eq. (|TTJ|) is 
the ith row of Eq. © and {y 5 }^! are independent, 
thus we can collect Eq. (fT0|) into i"C terms and the sth 
(s = 1, 2, • • • , if) term describes the degree of subsystem 
i contributed by sth cluster. If subsystem i belongs to 
cluster s', then the sth term will be 

N s +n s 

V s (Ci j -C mj ) = 0,(i = N s ,+l,N s ,+2,--- ,N s ,+n s ,) 

j=N„ + l 

(13) 

where N s — X^i=i n i an( i Tn = i + 1 when i < N s i + n s i ; 
m = 1 when i = _/V s / + n s / . 

Considering two different situations in Eq. (| 13[) . we 
can draw two statements as follow: 

SI: s ^ s'. ^Lijlpff+iCij is the degree of subsystem 
i, which is contributed by the sth cluster. So Eq. (fT3"|) 
shows that the degrees of subsystems in a cluster 
(e.g., the s'th cluster), which are contributed by the 
subsystems in another cluster (e.g., the sth cluster), 
should be the same. 

S2: s = s' . Since C mm — Cu = — 1, the rest part 
of Eq. (fT3"l) illustrates that the subsystems' degrees 
contributed by their cluster also should be the 
same. 

The two statements are the complete representation of 
Eq. ^ and now we have a clear physical picture of the 
necessary and sufficient condition of the existence of a 
PaS solution. Then, by using the two statements in the 
particular cases, we can obtain some interesting results: 

I. The nonexistence of any PaS states in the system 
shown in Fig. [T]is easy to be proved. Now we do not need 
to substitute all the possible solutions into the dynami- 
cal system. One can suppose that there is at least one 
PaS solution. According to SI and S2, the PaS should 
be observed between the subsystems with the same de- 
grees, so the subsystem 2 itself must be a PaS cluster 
without anyone else. On the other hand, subsystem 2 
is connected to subsystems 1, 3, 5, 6. SI requires that 
these four subsystems should be included in one or more 
clusters, but any of these clusters will not include sub- 
system 4, because the degree of subsystem 4 contributed 
by node 1 is while that of others is 1, i.e., subsystem 4 
itself should be a cluster too. And then, subsystem 4 is 
connected to subsystem 3 and 5, and because of the same 
reason, subsystem 3, 5 should be in different clusters with 
subsystem 1,6. However, the degrees of subsystem 3 and 
subsystem 1 are different with that of subsystem 5 and 
subsystem 6 respectively, that means they also should be 
in different clusters. So, we can see there will be at least 
6 clusters in this system with just 6 subsystems, any PaS 
phenomenon will not be observed in this system. 
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II. If there is symmetry in a dynamical system, Eq. 
([5]) can be rewritten as 

Cij — C mn , (14) 

where i = N a > + 1, N a > + 2, • • • , N a > + n s -\ j = N s + 1, 
N s + 2, ■ ■ ■ , N s + n s , N s = n i an d m = i + 1 when 

i < N s i + n s i ; m = 1 when i = iV s ; + n s / ; n = j + 1 when 
j < N s , + n s > ; n = 1 when j = iV s > + n s > . So Eq. (H]) is 
a stronger condition than Eq. (flU]) . That is why in few 
body systems @, H, there are tight relations between 
the PaS solution and the symmetry in the coupling struc- 
ture. Furthermore, one can prove that if there are only 
one or two subsystems in every cluster (n s i ,n s = 1 or 
2, s,s' = 1, 2, • • ■ , K) then there should be at least one 
symmetry in the coupling structure (n s i ,n s = 2 is non- 
trivial). In this case, the (Cij)rv xn s can only be ( ^ J J 

or { j J J and Eq. ([14")) will be satisfied automatically. 

III. The CS, as a special case of the PaS, can also be 
included in the Eq. ^ , where the dimension of V is only 
1. And its normalized basis vector is X 9 ; f, a i = ljy.i/V^V- 
In fact, this vector is always the eigvector of all of T. 
because every PaS subspace contains the CS subspace 
(x 1 = x 2 = • • • = x w ). The C is a non-row matrix (the 
summery of any row is 0) and the matrix product TC, CT 



also are. Thus, we conclude that, for the coupled system 
Eq. (JTJ) , the CS solution always exists no matter what the 
coupling structure C is. So Eq. ([9]) can be satisfied for 
all of C in this case, i.e., the CS solution exists in every 
system Eq. ([I]) no matter what the system topology C 
is. 

In conclusion, we have studied the relationship between 
the coupling structure and the PaS state in general dy- 
namical systems. The sufficient and necessary condition 
of the existence of PaS state for the coupling structure is 
found by the exact proof. And the result is counterintu- 
itive, that is, the existence of PaS state doesn't require 
the symmetry in the coupling structure, which is sup- 
posed before. Then, a new structure, the equal-degree 
random structure, is obtained. According to the suffi- 
cient and necessary condition, it is the general structure 
for the existence of the PaS state. Finally, we should 
stress that the proof also can be applied to the differen- 
tial dynamical systems and the conclusions are the same. 

The authors are grateful for discussions with Ping Lin 
and Jinfang Zhang. The work was in part supported 
by the 100 Person Project of the Chinese Academy of 
Sciences, the China National Natural Science Founda- 
tion with Grant No. 10775157, the Hong Kong Research 
Grants Council (RGC), and Hong Kong Baptist Univer- 
sity. 



[1] Y. Tanguy, et ai, Phys. Rev. Lett. 96, 053902(2006); A. 

Vilfan, F. Julicher, Phys. Rev. Lett. 96, 058102(2006); I. 

Fischer, et ai, Phys. Rev. Lett. 97, 123902(2006). 
[2] E. Emberly, N. S. Wingreen, Phys. Rev. Lett. 96, 

038303(2006); D. Battogtokh, K. Aihara, J. J. Tyson 

Phys. Rev. Lett. 96,148102(2006); Y-C. Lai, M. G. 

Frei, I. Osorio and L. Huang, Phys. Rev. Lett. 98, 

108102(2007). 

[3] R. E. Amritkar, G. Rangarajan, Phys. Rev. Lett. 96, 
258102(2006). 

[4] K. Kaneko, Phys. Rev. Lett. 65, 1391(1990); D. H. 
Zanette, A. S. Mikhailov, Phys. Rev. E 57, 276 (1998); 
O. Popovych, A. Pikovsky, and Y.Maistrenko, Physica 
D 168-169, 106 (2002); S. Sinha, Phys. Rev. E 66, 
016209 (2002); S. Jalan, R. E. Amritkar, GK. Hu, Phys. 
Rev. E 72, 016211 (2005); R. E. Amritkar, S. Jalan, 
GK. Hu, Phys. Rev. E 72, 016212 (2005). J. Garcia- 
Ojalvo, Y. Moreno, A. Arenas, Phys. Rev. Lett. 98, 
034101(2007); M. Rosenblum, A. Pikovsky, Phys. Rev. 
Lett. 98, 064101(2007); H. Nakao, K. Arai and Y. Kawa- 
mura, Phys. Rev. Lett. 98, 184101(2007). 

[5] This is one of widely used forms of the coupling non- 
linear maps. One can apply certain m-dimensional map 
x n +i = f(x„), x = (a;„,x^, ••■ ,x™) as subsystem to 
construct a coupled system X = (x\x 2 , • • ■ , x ), F = 
(f 1 ,f 2 ,--- , f ). The common interesting is the linear 
coupling between the subsytems, thus the coupling term 
will be eT ® CF(X), where e the coupling strength. V : 
R m — * R m characterizes the coupling scheme among the 



subsystems and C is the Laplacian matrix, dj = Aij/ki 
for j 7^ i and Cu — — 1, where ki is the degree of subsys- 
tem i and Aij is an element of the adjacent matrix A of 
the coupling structure. To gain insight in numerical simu- 
lation, we set f(x) to be the logistic map f(x) = 4a;(l— a;), 
but the proof do not depend on the dynamics. 

[6] L. M. Pecora and T. L. Carroll, Phys. Rev. Lett 80 
(1998) 2109; A. Pikovsky, M. Rosenblum and J. Kurths, 
Synchronization: A universal concept in nonlinear sci- 
ences (Cambridge University Press, Cambridge, 2001); 
C.S.Zhou et ai, Rev. Mod. Phys. (unpublished). 

[7] S. C. Manrubia and A. S. Mikhailov, Phys. Rev. E 60, 
1579(1999). 

[8] Y. Jiang et al, Phys. Rev. E 68, 065201 (R) (2003). 
[9] Y. Zhang et al, Phys. Rev. E 63, 026211(2001). 
[10] In this system, one can define the average dis- 
tance between two oscillators i and j by dij = 

limN^ooiN^ 1 J2n = l H X « _ X nll)> and SaV 8X6 in a 

synchronous state in the sense that dij < 10~ 12 . And the 
partial synchronization of a network means that some of 
nodes are in a synchronous state while others not. 

[11] R. Monasson, Eur. Phys. J. B 12, 555 (1999). 

[12] G. B. Dantzig, M.N.Thapa, Linear Programming P.327, 
(Springer- Verlag, 1997); L. Smith, Linear Algebra P. 156, 
(Springer, 1998). 

[13] The Td in the second example is only a special case of 
Eq. (|12p . and the permutation transformation done to A4 
can also be represented by T4 = R2 © 1 ffi 1. 



